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ABSTRACT 

In this paper we report on the early evolution of a core-collapse supernova explosion 
following the birth of a magnetar with the dipolar magnetic field of _B = and 
the rotational period of 2ms, which was studied by means of axisymmetric general 
relativistic MHD simulations. In this study we use realistic EOS and take into ac- 
count the cooling and heating associated with emission, absorption, annihilation, and 
scattering of neutrinos, the neutrino transport is treated in the optically-thin regime. 
The supernova explosion is initiated via introducing into the initial solution the "ra- 
diation bubble" , whose total thermal energy is comparable with the typical energy of 
supernova ejecta. The numerical models exhibit highly coUimated magnetically-driven 
jets very early on. The jets are super-Alfvenic but remain sub-fast until the end of 
the simulations (t=0.2s). The power released in the jets is about 3 x 10^°erg/s which 
implies the spin-down time of ~ 37s. The total rotational energy of the magnetar, 
E ~ lO^^erg, is sufficient to drive a hypernova but it is not clear as to how large 
a fraction of this energy can be transfered to the stellar ejecta. Given the observed 
propagation speed of the jets, Vp ~ 0.17c, they are expected to traverse the progenitor 
in few seconds and after this most of the released rotational energy would be sim- 
ply carried away by these jets into the surrounding space. 3-dimensional effects such 
as the kink mode instability may reduce the jet propagation speed and increase the 
amount of energy transferred by the jets to the supernova ejecta. Our results provide 
the first more or less self-consistent numerical model of a central engine capable of 
producing, in the supernova setting and on a long-term basis, coUimated jets with suf- 
ficient power to explain long duration GRBs and their afterglows. Although the flow 
speed of our jets is relatively low, only Vj ~ 0.5c, the cooling of proto-neutron star 
will eventually result in much higher magnetization of its magnetosphere and ultra- 
relativistic asymptotic speeds of the jets. Given the relatively long cooling time-scale 
we still expect the jets to be only weakly relativistic by the time of break out. This 
leads to a model of GRB jets with systematic longitudinal variation of Lorentz factor 
which may have specific observational signatures both in the prompt and the afterglow 
emission. The simulations also reveal quasi-periodic ejection of plasma clouds into the 
jet on a time-scale of 20ms related to the large-scale global oscillation of magnetar's 
magnetosphere caused by the opening-closing of the dead zone field lines. These kind 
of central engine variability may be partly responsible for the internal shocks of GRB 
jets and the short-time variability of their gamma-ray emission. 

Key words: supernovae: general - stars: neutron - gamma-rays: bursts - methods: 
numerical - MHD - relativity 



1 INTRODUCTION 

After decades of intensive research the exact mechanism of 
core-collapse supernovae (SNe) still remains a mystery. It 
is widely believed that the explosion is driven by neutrinos 
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(SSK); 



emitted by the proto-neutron star (PNS) formed as the re- 



sult of the collapse (Bethe 19901. However, the attempts to 



reproduce core-collapse SNe in computer simulations have 
encountered severe problems. The reason for the failures be- 
comes more or less clear when one compares the total energy 
radiated in the form of neutrinos, ~ 10^'^erg, with the typ- 
ical energy of supernovae, ~ lO^^erg. This tells us that in 
order to obtain a reliable answer the computational error in 
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the energy budget has to be below 1%, which a is very de- 
manding condition. The neutrino transport, of six different 
species, has to be treated with great accuracy including the 
processes effecting heating and cooling rates. 

The alternative magnetic mechanism of core-collapse 
SNe has been around for a while, first proposed by 
Bisnovatyi-Kogan(p970| and LeBlanc & WilsonJ 19701), 



whose numerical simulations were miles ahead of their time. 
For three decades this mechanism was not taken seriously 
and only occasionally efforts were made to develop it fur- 



of the black hole and/or the accretion disk (e.g. Lyutikov & 
Blandford 2003, Vlahakis & Konigl 2003, Proga et al., 2003). 

The energy of HNe ejecta is of the same magnitude 
as the rotational energy of a neutron star (NS) with very 
short period of rotation, 1-2 milliseconds. This energy can 
be released via magnetic braking within a very short period 
of time, from few seconds to few hundreds of a second, pro- 
vided the magnetic field is very strong, of order _Bns — lO'^^G 



ther ( Bisnovatyi-Kogan 1976 Meier et al. 1976 Symbalisty 



1984|). Now it is experiencing "renaissance" (Wheeler et al 
2000; Wheeler et al. 2002[ [Akiyama et al. 20031 Yanmda'fc 



Sawai 2004, ,Kotake et al. 2004| |Mizuno et al. 2004, ^Taki- 



2005 



2006 



2007 



waki et al. 2004| Xrdeljan et al. 2005| Akiyama fc Wheeler 



Proga 20051 'Moiseenko et al. 2006 [Obergaulinger 



Shibataet al. 2006 , Masada et al. 2007| |Nagataki ct_al 



( jThompson 2006 Bucciantini et al. 2006 1. This value is three 
orders of magnitude above the typical strength of magnetic 
field of normal neutron stars and would be considered a huge 
over-stretch not so long time ago. However, the identifica- 
tion of Soft Gamma Ray Repeaters and X-ray Anomalous 
Pulsars with magnet ars, neutron stars with magnetic field 



10'" - lO'^^G (Thompson & Duncan 1995 



Woods & 



Thompson 20041, shows that the magnetar model may be a 



viable alternative to the coUapsar model of HNe as well as 



Burrows et al. 20071, thanks to the slow and difficult 



progress of the neutrino models of SNe, the development of 
robust numerical methods for MHD, the accumulated evi- 



GRBs (Usov 1992 Thompson et al. 2004 Thompson 2006 



Metzger et al. 2007 Uzdensky fc MacFadyen 2006 1 



deuce for asphericity of supernovae (Khokhlov et al. 1999 
Wang et al. 2002[[Wang et al. 2003| P] the high coUimatiori 



of flows associated with GRBs (e.g. Piran 2005a,b) and the 
increasing popularity of the magnetic mechanism for the ori- 
gin of other astrophysical jets. The results show that mag- 
netic field can facilitate explosions of standard power and 
even more powerful explosions provided the magnetic field 
is sufficiently strong. 

The energy of high-velocity supernovae or hypernovae 
(HNe) exceeds that of normal core-collapse supernovae by 
about ten times (e.g. Nomoto et al. 2004). The few SNe 
that have been reliably identified with GRBs are all HNe. 
Although, the statistics is very poor and is subject to ob- 
servational bias towards most powerful and hence brightest 
events, the connection between HNe and GRBs is widely ac- 
cepted. The physical mechanism of HNe is not established 
and there is no shortage of competing theoretical models. 
The most popular one, the collapsar model, relates HNe 



The origin of magnetar's magnetic field is still debated. 
One possibility is the generation of unusually strong mag- 
netic field in convective cores of some pre-supernova stars. 
During the core collapse the poloidal field is amplified like 



and in order to reach the required Bn 



lO'^'G af- 



ter the collapse the original field in the core, on the scale 
of lO^km, should be of order Be ~ 10"G. The coUapse of 
magnetized rotating stellar cores has been studied by many 
groups. The general conclusion from the latest parameter 
studies is that only for Be > lO'^G the magnetic field is 
strong enough to influence the dynamics of supernovae ex- 
plosions ( [Obergaulin ger 2006! [Burrows et al. 2007[ ). The key 
feature of these explosions is strong differential rotation and 
generation of toroidal magnetic held whose pressure actually 
drives the explosion. The hoop stress of this field also en- 
sures anisotropy of the explosion - it is more powerful along 
the rotational axis. 

On the other hand, the super-strong magnetic field of 
magnetars could be generated via Q-Sl-dynamo in the con- 



with collapse of m assive rapidly rotating stars ( |MacFadyen| vective PN S puncan fc Thompson 1992[ [Thompson fc Dun 
fc Woosley 1999[ ). In this model, the stellar core collapses [can 1993[ ). In this theory the strength of saturated field 

strongly depends on the rotational period - shorter period 
leading to stronger field. In order to generate the dipolar 
field of Bns ^ 10"G the period of PNS should be very 



straight into a black hole (which used to be considered as 
an indication of "failed" supernova explosion) but the stellar 
envelope does not and forms a hyper-accreting disk around 
the hole. The energy released in the disk can be very large, 
sufficient to drive HN and GRB, but the way it is trans- 
ferred to the supernova ejecta and the GRB jets remains 
to be determined. It is widely assumed that GRB jets are 
powered by the energy released in the funnel region of the 



short, not much higher than few milliseconds (Duncan fc 



Thompson 19921. Thus, this theory unites both basic con- 



disk via annihilation of neutrinos (Ruffert fc Janka 1999 



Narayan et al. 20011. However, taking into account opacity 



in the inner regions of the disk significantly reduces the ef- 



ficiency of this mechanism ( Di Matteo et al. 2002 Nagataki 



et al. 20071. Alternatively, as it has been suggested by many 



authors, the GRB-jets in the collapsar scenarion can be pow- 
ered via magnetic mechanism tapping the rotational energy 



^ Based on the results of 2D axisymmetric numerical simulations, 
it has been proposed that strongly aspherical explosions may arise 
in non-magnetic models too (Blondin et al. 2003 Burrows et al. 



ditions required to produce a hypernova in the magnetar 
model, rapid rotation and super-strong magnetic field, in 
one. On the downside, the PNS is no longer magnetically 
coupled to the envelope and the transfer of its rotational 
energy to the stellar envelope is more problematic. Turbu- 
lence required for the magnetic dynamo action can also be 
generated via the magneto-rotational instability (MRI, Bal- 
bus fc Hawley,1991). Calculations based on the linear theory 
show that in the supernova context strong saturation field 
can be reached very quickly on the time scale of only several 



20061. 3D simulations are needed to test this idea. 



tens of rotational periods (Akiyama et al. 20031. 

Thompson et al.(2004) proposed a model of magneti- 
cally driven HN explosions and GRB where the usual de- 
layed neutrino-driven mechanism plays the role of bomb 
detonator. The neutrino-driven blast creates a rarefaction 
around the newly-born magnetar, its magnetosphere ex- 
pands and eventually develops a magnetically driven pulsar 
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wind (PW). This wind acts as a piston on the expanding 
supernova envelope, that is the layer of compressed stellar 
material behind the supernova shock, and energizes it be- 
yond the level of normal SN. 



Metzger et al. 1 2007 1 generalized the ID model of equa- 
torial magnetically-driven wind by Weber & Davis ( [T967[ ) to 
fit the conditions of PNS. Axisymmetric magnetized winds 
of magnetars with various degree of mass loading have been 
studied in Bucciantini et al. (2006J via 2D numerical simula- 
tions. They concluded that during the initial phase (Kelvin- 
Helmholtz timescale) when the wind is already magnetically- 
driven but still non-relativistic the rotational energy loss 
rate is L ~ 4 x 10^^-Bi5Pj~^^^erg/s, where -B15 is the mag- 
netic field of the PNS in lO^^G and Pi is its period in mil- 
liseconds. Bucciantini et al.( |2007a[ ) applied the model of pul- 
sar wind nebula (PWN) by Begelman & Li ( |1992[ ) to study 
the late phases of hypernova explosion, t > Is. They used 
the thin-shell approximation to describe the dynamics of 
the supernova shock driven by the anisotropic pressure of 
the magnetized bubble, "baby pulsar wind nebula" , created 
inside the exploding star by the magnetar wind. They con- 
cluded that when the magnetic energy in the bubble exceeds 
~ 20% of the thermal energy the shell expansion proceeds in 
a highly anisotropic fashion and eventually an axial channel 
is created inside the stellar envelope. They proposed that 
later this channel coUimates the ultra-relativistic wind from 
cooled magnetars into a pair of GRB jets. In the follow up 



paper ( Bucciantini et al. 2007b I they studied the interac- 



tion of the super-fast magnetar wind with the supernova 
envelope via fully relativistic MHD simulations and arrived 
to similar conclusions. For the wind model to become ap- 
plicable the size of the "central cavity" must exceed, and 
likely many times, that of the fast surface as found in the 
unconfined case. Prior to this one would expect a somewhat 
different dynamics. 

The most pessimistic possibility is a rigidly rotating 
magnetosphere in which case no extraction of rotational en- 
ergy takes place. However, the theory of relativity places a 
upper limit on the size of such a magnetosphere. Indeed, the 
magnetosphere whose size exceeds the light cylinder (LC) ra- 
dius cannot be rigidly rotating as this would imply superlu- 
minal motion of magnetospheric plasma and should become 
differentially rotating. This limit is particularly relevant for 
the magnetospheres of millisecond pulsars as the radius of 
LC, TIlc — 50Pi km, is only few times the radius of the 
star itself. Once the differential rotation sets up it amplifies 
the toroidal field, the magnetic pressure grows and promotes 
further expansion of the magnetosphere. Uzdensky & Mac- 
Fadyen([2006) suggested that this factor along can explain 
supernova explosions but they did not elaborate ways to 
achieve the initial expansion beyond LC. Moreover, they as- 
sumed low mass loading of magnetic field lines, more suitable 
for NS than hot proto-NS, and concluded that the magneto- 
spheric plasma will exhibit ultra-relativistic rotation in the 
equatorial plane. In any case, this argument suggests that 
the magnetic braking should begin to operate at the very 
early stages of the "detonated" magnetically-driven explo- 
sions, well before the central cavity expands beyond the fast 
surface of unconfined magnetar wind and this may to have 
important implications for their later development. 

In this paper we describe an attempt to explore the 
early phase {t < Is) of such explosions via axisymmetric 



General Relativistic MHD simulations. Our main goals are 
to study the basic dynamics of this phase, to determine the 
spin-down power, and to assess the implications for the later 
evolution. We use an upwind conservative scheme that based 
on a linear Riemann solver and uses the constrained trans- 
port method to evolve the magnetic field. The details of 
this numerical method and various tests are described in 
Komissarov(!l999'; '2004b: 2006). In Sec.2 we describe the de- 
tails of numerical experiments, including initial and bound- 
ary conditions. Sec. 3 describes the results of simulations, 
and in Sec. 4 we discuss their implications for the model of 
magnetar-driven supernovae explosions and GRB jets. Our 
conclusions are listed in Sec. 5. 



2 SIMULATION SETUP 

The ultimate hypernovae/supernovae simulations should 
trace all phases the stellar explosions, the onset of core col- 
lapse, core bounce, delayed explosion, generation of mag- 
netic field, formation of magnetic cavity, and break out of 
magnetically-driven jets. Unfortunately this is still beyond 
our current capabilities and we have to make some simplifi- 
cations. Here we start our simulations from the point when 
strong magnetic field has just appeared from surface of al- 
ready formed PNS. To account for the previous stages of 
the explosion we assume that PNS is surrounded by a "ra- 
diation bubble" with a large amount of energy deposited by 
neutrino in form of heat. The bubble is surrounded by the 
collapsing stellar envelope. 

The gravitational attraction of PNS is introduced 
via Schwarzschild metric in Boyer-Lindquist coordinates, 
{(ji,r,0}. The two-dimensional computational domain is 
(ro < r < ri) X {0 < 6 < tt), where tq corresponds to 
the radius of PNS and ri — lO^km. The total mass within 
the domain is small compared to the mass of PNS that al- 
lows us to ignore its self-gravity. The grid is uniform in 6 
where it has 200 cells and almost uniform in log(r) where 
it has 430 cells, the linear cell size being the same in both 
directions. 



2.1 PNS 

The evolution of a non-rotating proto-neutron star of mass 
M — 1.4 M0 h ave been studied in details by Burrows & 
Lattimer (19861. According to this study the deleptoniza- 



tion time-scale is rather long, it takes around 5 seconds be- 
fore half the leptons inside M{r) < 0.5M are lost. However, 
the contraction time-scale is much shorter - it takes only 
~ 250ms for the star radius to decrease down to ~ 15km 
(see fig. 5 in Burrows & Lattimer, 1986). Using these re- 
sults we fix the mass and radius of PNS to M = 1.4Mq 
and ro — 15km respectively, and we choose the rotation pe- 
riod to be P = 2ms. To set the boundary conditions on 
the PNS surface we utilize the results of high resolution ID 



numerical models of PNS winds ( Thompson et al. 2001 Met- 
|zger et al. 2007[ ). The gas temperature of "ghost cells" is to 
To = 4 MeV, which is typical for a hot newly-born PNS, 
and their density is po = 3 x lO^gcm"'^. This density cor- 
responds to the interface between the very thin exponential 
atmosphere of PNS and its wind. The initial magnetic field 
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is purely poloidal with the azimuthal component of vector 
potential 



where 



fir) 



l-((r-ro)/(re-r„))2 if 
if 



r < rc 
r > rc, 



(1) 



(2) 



Bo = lO^^G and = l.SrcQ For /(r) = 1 equation [T] de- 
scribes magnetic dipole of strength Bo at r = ro, S = 0. 
The masking function ^ "moves" the magnetic field lines 
inside the sphere of radius Vc without changing the mag- 
netic flux distribution over the PNS surface. This distribu- 
tion is preserved for the whole duration of simulations. The 
radial velocity of ghost cells is set to zero and the mass out- 
flow through the boundary is found via solving the Riemann 
problem at the boundary interface. 

We assume that neutrino luminosity in each species is 
given by the black body formula 



.20 X W^\rlc ( -^Y ^ ^ 6.5 X 10' 
° VlMeVy s 



51 erg 
s 



(3) 



and that the mean neutrino energy is = 3.15 To — 
12.6 MeV. This oversimplification is unlikely to have a strong 
effect on the largely magnetically driven outflows from PNS. 



2.2 Collapsing star 



Here we adopt the simple free-fall model by Bethe (19901. 
According to this model the radial velocity is 

= (2GM/^)^/^ 

and the density is 

— '°'(ll)"'(lo^)"'^ 

where Ci is a coefficient between 1 and 10 ( [Bethe 1990[ ) and 
ts is the time since the onset of collapse. The corresponding 
accretion rate and ram pressure are 

1/2 



(4) 



(5) 



M = 0.038 Ci 



M 



I AM, 



3.7x10^*^(71 



I AM, 



^ (r^) (lO^cm) 



5/2 



(6) 



r(7) 



respectively. Since in the simulations we flx the core mass to 
M = 1.4Mq the model is fully determined by the combina- 
tion C\/ts. For a core of radius rc = lO^cm and mass 1.4A'^q 
the collapse duration can be estimated as tc — 2rc/3i;g- ~ Is. 
The time-scale of the delayed explosion is likely to be around 



few tens of a second after bounce (Bethe 19901. Thus, the 
value of ts is relatively well constrained and we fixed it to 
be ts — Is. Since GRBs are currently associated with more 
massive progenitors we consider Ci = 3,9. 

To account for rotation of the pre-supernova star the 
accreting mass is attributed with specific angular 



^ Here we give the component of vector potential in non- 
normalized coordinate basis. Everywhere else in the paper the 
components of vectors are measured using normalized bases. 



Model 


A 


B 


C 


D 


Ci 


3 


3 


9 


9 




1 


0.1 


1 


0.1 



Table 1. Parameters of the models used for test simulations. 



I = lo sin 



(8) 



where lo — lO^^'cm^s^^. This is slightly higher then the 
equatorial value of I on the PNS surface, ~ 7 x lO^^cm'^s"^ 
but still low enough to prevent the development of accre- 
tion disk around the magnetar. In fact, the rotation does 
not seem to have a strong effect on the solution. 

Eqs.( 4|5p l are also used to set the flow variables in the 
ghost cells of the outer boundary. 



2.3 Radiation bubble 

The radiation bubble is assumed to extend up to rt = 
300 km or ri, = 200 km in different models. The radial and 
polar velocity components of matter in the bubble are set to 
zero and the azimuthal component is calculated using eq. (|8| . 
The prescribed density and pressure distributions are 



Pb 



(^) 



^ro 
GMpt 



-(l + a) 



1 



+ 



2A GM 



l + a 2 + a rc^ 



Pb- 



(9) 
(10) 



For A = 1 these distributions correspond to a hydrostatic 
equilibrium. In the simulations we use pt = O.OSpo and vary 
Pb to get the desired amount of thermal energy stored in 
the bubble, Ei, = CalO^^erg. The parameters of four models 
described in this paper are summarized in Table 1. 



2.4 Microphysics 

We use realistic equation of state (EOS) that takes 
into account the contributions from radiation, lep- 
ton gas including pair plasma, and non-degenerate 
nuclei - we have incorporated the EOS code HELM 
kindly provided by Frank Timmes for free download 
( http://www.cococubed.com/code-pages/eos.shtmll . 
In this code the contribution from the electron-positron 
plasma is computed via table interpolation of the Helmholtz 
free energy ( [Timmes fc Swesty 2000"[ ). The neutrino trans- 
port is treated in the optically thin regime. The neutrino 
cooling is computed using the interpolation formulas given 
in Ivanova et al.( [Ivanova et al. 1969[ ), for the Urca-process, 
and in Schinder at al.(1987l, for pair annihilation , photo- 
production, and plasma emission. The neutrino heating 



rates are computed using the prescriptions of Thompson 



et al.(2001| which take into account gravitational redshift 
and geodesic bending. We ignore the Doppler effect due 
to plasma motion as its speed relative to the grid never 
becomes highly relativistic (see Sec|3|. Photo-disintegration 
of nuclei is included via modifying EOS (e.g. Ardeljan et 
al.2005). The equation for mass fraction of free nucleons is 



adopted from Woosley & Baron ( 1992 1 
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Figure 1. Model A at time t ~ 200 ms. Left panel: logiop measured in g cm^''. Right panel: Radial velocity, v'~ /c. The unit length in 
all figures in this paper is £ = GM/c^ ~ 2km. The dynamic range of colour plots does not always reflect the full range of variation of 
the represented quantity but is rather selected to make more revealing images. 
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Figure 2. Same as in figure [l]but for model B at i ~ 200 ms. 
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3 RESULTS 



Figure [T] shows the global structure of model A at time 
t ~ 200 ms when the simulations were terminated because 
the shock wave has reached the outer boundary of the com- 
putational domain. The low density "column" along the 
symmetry axis reveals the volume occupied by the coUi- 
mated outflows from the magnetar. One can see that these 
jets have already "drilled" holes in the supernova envelope, 
that is the layer of stellar material compressed by the super- 



nova shock, and are beginning to propagate directly through 
the collapsing star. The unit length in this and other figures 
is C = GM/(? ~ 2km. Thus, the propagation speed of the 
jets is about 5 x lO^kms"'^ or 0.17c. At this speed the jets 
would travel across the star of radius ~ 2 x lO'^km in about 
4 s. The right panel of fig jl] shows that the flow speed of the 
jets can be as high as 0.5c. The plots show no signs of the 
jet termination shock - this is because the jet flow is sub- fast 
and can decelerate smoothly when it reaches the jet head. 
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Figure 3. Same as in figures^ 

The knotty structure of jet is not due to internal conical 
shocks or the pinch instability but due to the non-stationary 
nature of the central engine. The dead zone of magnetar's 
magnetosphere opens up in a quasi-periodic fashion with a 
characteristic time of 20ms. This causes significant restruc- 
turing of the open magnetosphere and variability of the out- 
flow. This is also the origin of waves clearly visible in the 
velocity field within the circular region of 3000km radius 
(right panel of figjl]) . One can also see weaker outflow along 
the equatorial current sheet of magnetars magnetosphere. 

Figure [2] shows the global structure of model B at the 
same time. The most apparent difference is in the size of 
the quasi-spherical supernova shock. It propagates notice- 
able slower because of the lower energy deposited in the ra- 
diation bubble. However, the propagation speed of the jets is 
almost the same as in model A suggesting that their power 
must be similar as well. One can also see that inside the su- 
pernova envelope the jets propagate inside a relatively thin 
channel and outside of the envelope they propagate through 
a low density cocoon inflated during the previous evolution. 
A similar cocoon is beginning to form in model A (see flgjl]) 
suggesting that at later times models A and B could look 
even more alike. 

Figure [3] shows the global structure of model C at time 
t = 277ms. Since the ram pressure of the infalling material is 
higher than in models A and B the jet propagation speed is 
somewhat lower. The same factor explains lower propagation 
speed of the supernova shock in this model compared to 
model A. In other respects the solution is not much different. 

The evolution of model D starts in very much the same 
fashion as for other models but at later times we observe 
a substantial fall-back of the shocked stellar material on 
the magnetar. The neutrino-cooled accreting matter even- 
tually blocks the jet channel and squashes the magneto- 
sphere against the star surface (see flgEl. Soon after the 
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but for model C at t ~ 277 ms. 




-10 10 



V B 9 10 11 

Figure 4. Model D at time t = 169ms. The colour image shows 
logiop and the contours show the magnetic field lines. The unit 
length in all figures in this paper is £ = GMjc? ~ 2km. 



code crashes as the conditions near the surface become too 
extreme. 

Figure [5] shows the total flux of "free energy", the total 
energy at infinity minus the rest mass energy, as a function 
of radius for models A,B, and C. The "spikes" at r ~ 1000 
for model B, between r = 1000 and r = 4000 for model 
A, and between r = 1000 and r = 3000 for model C cor- 
respond to the blast wave of normal supernova explosion. 
The "plateaus" to the left and to the right of the spikes 
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Figure 5. Total free energy flux through sphere of radius r for 
model A (solid line), model B (dashed line), and model C (dash- 
dotted line). 



correspond to the jets. Thus, in all models the spin-down 
power released in the jets is indeed more or less the same, 
L ~ 3 X lO^^erg/s. The corresponding spin-down time is 
~ 37s suggesting that most of the rotational energy will be 
released after the jets break out from the progenitor star. 
While this is a promising result for the magnetar model of 
GRBs it also suggests that only a modest fraction of the total 
rotational energy of the magnetar can be transferred to the 
supernova ejecta. In all three models the ratio of Poynting 
flux to the rest mass energy flux is a ~ 2.5. This is ap- 
proximately 5 times higher than the value predicted by the 
theory of unconfined equatorial magnetar wind (see eq.(27) 
in Metzger et al.(2007|). 

Figure [6] shows model A at the early times, t ~ 21 ms. 
One can see that the magnetar is already producing an out- 
flow which has blown a cavity in the shape of a column in- 
side the collapsing star. The radius of this column is close to 
the radius of magnetar's Hght cylinder, _Rlc — SOMG/c^ = 
100 km, which is shown in the right panel of fig[6] by the 
dashed line. We note, however, that for z < 100 the col- 
umn is still inside the light cylinder. The total pressure of 
the outflow is dominated by the magnetic field, the ratio of 
magnetic to gas pressure varying between 10 and 300. One 
can see strong axial compression of the outfiow due to the 
hoop stress of azimuthal magnetic field. The distribution of 
total pressure over the boundary of the volume occupied 
by the outflow has a clear maximum at the point of in- 
tersection with the polar axis which is obviously related to 
the axial compression. This distribution explains the faster 
expansion of the cavity blown by the outflow in the axial 
direction. In the middle panel of flg|6] one can see the bow 
shock, at 2 ~ 600, driven by the onflow into the surround- 
ing. There are no indications of the termination shock of the 
outflow itself - this is because the outflow is sub-fast. The 
solid line in the right panel of flg|6] shows the locus of points 
where the radial component of the phase velocity of ingo- 



ing Alfven waves becomes positive. Inside the column these 
point form a closed surface which we call the Alfven hori- 
zon of the magnetar-driven outflow. Any Alfven wave can 
propagate from the exterior of this surface into its interior. 
One can see that the Alfven horizon is always located inside 
the light cylinder, which agrees with theory of steady state 
MHD flows ( [Okamoto 1978 1. 

Figure [T] shows the inner region of model A by the end 
of the run {t ~ 200 ms) . By this time the size of the central 
cavity, that develops a diamond-like shape, is significantly 
larger than -Rlc, thanks to the decreased pressure in the 
surrounding bubble that resulted from its expansion in the 
course of supernova evolution. However, the rotational veloc- 
ity in the equatorial plane is not ultra-relativistic indicating 
that the relativistic winding mechanism is not in operation. 
The velocity field (left and right panels of fig[7| shows that 
close to the magnetar the outflow is similar to an isotropic 
wind but then it becomes progressively coUimated in the 
polar direction and eventually passes through the "bottle- 
neck" nozzles at the axial corners of the diamond. The coUi- 
mation is smooth with no indication of shocks waves which 
is expected because the flow is always sub-fast. 

After passing trough the nozzle each jet enters the chan- 
nel made during the previous evolution and propagates in- 
side a low density cocoon. In its outer layer the cocoon con- 
tains poloidal field lines of opposite direction to that of the 
jet and closer to the jet one can see magnetic islands and ev- 
idence of vortex motion. These islands are mostly remnants 
of plasmons injected from the dead zone of magnetar's mag- 
netosphere (see the right panel of fig[8| during its oscillatory 
cycle. In the middle panel of flg[7] one can see one of such 
recently ejected plasmons. Close inspection of animated pic- 
tures suggests the following origin of the oscillations. When 
the dead zone includes the highest magnetic flux it grad- 
ually accumulates plasma supplied from the magnetar sur- 
face. This leads to increase of the centrifugal force and ex- 
panding of the dead zone. Its magnetic fleld lines stretch 
in the equatorial direction trying to open up. The process 
seems to accelerate as the dead zone plasma moves further 
out from the star and flnally escapes in eruptive manner 
leaving behind thin equatorial current sheet. Then the op- 
positely directed fleld lines of the sheet reconnect and this 
restores the initial configuration of the dead zone (a similar 
behavior of PNS magnetosphere was reported in Bucciantini 
et al., 2006). The amplitude of these oscillations is large 
so one gets the impression of the magnetosphere "breath- 
ing heavily" . Since here we integrate the equations of ideal 
RMHD it is the numerical resistivity that is responsible for 
the observed reconnection and the reliability of its time-scale 
is unclear. 

The plots of flg[7] also reveal a rather peculiar time- 
dependent structure inside the diamond cavity and near to 
the polar axis, to which we will refer as the "trap zone". It 
begins at 2 < 50 and continues almost up to the top corner of 
the cavity. In the right panel of fig[7[ that shows the ration 
of magnetic to gas pressure (the inverse of the traditional 
magnetization parameter 13), it looks as the region of lower 
magnetization. The density and velocity plots of fig[7]show 
that the trap zone has a relatively high density and low 
speed. In fact, the radial velocity in the trap zone changes 
sign both in space and time - the density clumps seen in the 
trap zone exhibit oscillatory motion along the polar axis. 
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Figure 6. Model A at 4 ~ 21 ms. Left panel: logj^Q p and the poloidal velocity field; Middle panel: The colour image shows log^gPtoti total 
pressure, and the contours show the magnetic flux function; Right panel: The colour image shows logjQ(p„/pg), the ratio of magnetic 
pressure to the gas pressure. The solid lines show the surfaces where the radial phase speed of ingoing Alfven waves is zero - the one 
closest to the symmetry axis is the "Alfven horizon" . The dashed line shows the location of the light cylinder. 



This explains tiie peculiar turns of magnetic field lines in 
the trap zone (see the middle panel of fig[7| . From time to 
time plasma clouds are ejected from this zone and the time- 
scale of this ejections is similar to the time-scale of global 
magnetospheric oscillations. One of such clouds is seen in 
density plot of fig[7]at z ~ 650. Another peculiarity of this 
zone, to which we have no explanation, is that it is separated 
from the polar axis by cylindrical shell, of radius R ~ 7?lc, 
with a relatively high velocity. 



Figure [S] provides with additional information on the 
magnetic structure of this solution. One can see that the 
azimuthal magnetic field dominates almost everywhere in 
the "diamond zone" and in the jet, with the exception of 
the region immediately around the symmetry axis. /Bp is 
particularly high in the middle of the magnetic cocoon where 
the poloidal field changes direction (see the left panel of 
fig|8|. However, the azimuthal field changes sign only in the 
equatorial current sheet. This is clearly seen in the middle 
panel of fig [sjwhich shows that B"^ does not vanish anywhere 
else. The electric current is concentrated in the jet core, 
the cocoon surface, and the equatorial current sheet thus 
showing that within most of the jet volume the magnetic 
field is almost force-free. The right panel of fig[8] shows the 
most inner part of the magnetar magnetosphere and reveals 
the existence of its dead zone. 



4 DISCUSSION 

4.1 Cautionary notes 

The main motivation of this study was to explore the effects 
that the birth of a millisecond magnetar may have on the 
development of delayed supernova explosions. The problem 
of supernova explosions is extremely complicated and we 
had to make a number of simplifying assumptions driven by 
the limitations of our numerical method. The most prob- 
lematic issue is the initial setup where we assume that the 
magnetic field of the magnetar generated via a-fi-dynamo 
or other dynamo has just emerged above magnetar's surface 
and the large amount of energy needed to drive standard 
supernova explosion have just been deposited in the radia- 
tion bubble. Although the time-scales of these processes are 
similar (< Is), even small asynchronism can significantly 
change the initial conditions simply because of the large ex- 
pansion speed of the supernova shock. This has to be kept 
in mind when analyzing our numerical models, particularly 
their early evolution. However, we expect the models to be 
more or less reliable in the second half of the simulation 
runs. 

Another important limitation of our models, as well 
as many other models of magnetically-driven astrophysi- 
cal fiows, concerns the potentially destructive role of the 
kink mode instabilities. Obviously, these instabilities are to- 
tally suppressed in our axisymmetric simulations and full 
3D simulations are needed to investigate this issue. We can 
only comment on two stabilizing factors present in the cur- 
rent setting. Close to the magnetar the poloidal magnetic 
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Figure 7. Model A at t ~ 200 ms. Left panel: logj^Q p and the poloidal velocity field; Middle panel: The colour image shows log^gPtot, 
total pressure, and the contours show the magnetic flux function; Right panel: log]^Q(p„/pg), the ratio of magnetic pressure to the gas 
pressure and the velocity field. 
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Figure 8. Model A at t ~ 200 ms. Left panel: The colour image shows logigQB'^]/ Bp^i), the ratio of azimuthal and poloidal components 
of magnetic field, and the contours show the lines of poloidal magnetic field; Middle panel: The colour image shows the magnitude of the 
azimuthal magnetic field, logj^Q and the contours show the lines of electric current. Right panel: The colour image shows the amplitude 

of B'f' and the contours show the lines of poloidal magnetic field. The unit length in all figures in this paper is C = GMjc? ~ 2km. 



field may provide strong "backbone" support for tlie out- 
flows (Anderson et al. 20061). Further away the jets propa- 



gate within a channel with "high-inertia walls" which could 
effectively dump the kink-type motions inside the channel. 

Finally, the condition of axisymmetry means that we 
can only study the exceptional case of aligned rotator. How- 
ever, we do not expect the case of oblique rotator to give 
dramatically different results. For example, in the limit of 
Magnetodynamics (inertia-free Relativistic MHD) the spin- 



down power of orthogonal rotator exceeds that of aligned 



rotator but only by a factor of two ( Spitkovsky 2006 1 



4.2 Comments on the mechanism of explosion 

Many features of the magnetar driven outflows observed in 
our simulations are similar to those predicted by Uzdensky 



& McFadyen (20061. There is however a number of impor- 



tant differences and one of them is the nature of explosion. 
None of the succesful explosions in our simulations are mag- 
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netically driven. These explosions are produced via the de- 
position of thermal energy in the radiation bubble of the 
initial solution. We have never seen the long phase (many 
rotational periods) of building up of magnetic pressure be- 
hind the stalled shock - the key process in their scenario. 
Moreover, in the simulation runs where the amount of en- 
ergy deposited in the buble was reduced in order to ap- 
proach the conditions of failed supernova and to engage the 
magnetic mechanism the magnetar magnetosphere did not 
expand but was pressed back to the magnetar surface un- 
der the weight of neutrino-cooled layers of accreting mate- 
rial. It is easy to see why. Suppose that mass AM passes 
through the stalled shock of failed supernova and accumu- 
lates at the radius of the light cylinder Tlc ~ 100km (the 
critical radius in the magnetic mechanism). Then its weight 
will be iV = GMAM/Anrtc =^ 3 x lO^^gcm'^s"^ for the 
relatively small AM = O.OIMq. This is much larger than 
the pressure of the dipolar magnetic field at the same ra- 
dius, Pm — 1.3 X 10^®gcm~^s~^ (for this estimate we as- 
sume that the at the PNS surface Bo ~ lO^'^G). Similarly, 
one finds that the ram pressure of the free-falling stellar 
material at the light cylinder is also much larger than the 
magnetic pressure. In principle, the magnetic field generated 
in the PNS can be transported in the region bounded by the 
stalled shock of failed supernova by the large scale convec- 
tive motions developing in this region and even to be locally 
amplified via some sort of dynamo. However, in this case the 
dynamics of magnetic field would be much more complicated 
than the one described in the simple magnetospheric model 
of Uzdensky & McFadyen {2006} and in order to capture 
it one would require much more sophysticated simulations 
than ours. 



4.3 Magnetised winds and magnetic towers 

Until recently the magnetically driven outfiows from astro- 
physical objects were classified only as "winds" or "jets", 
depending on whether the outfiow proceeds in all direc- 
tions or it is coUimated within a cone due external con- 
finement of one sort or another. But since the work of 



Lynden-Bell( 1996 20031, who proposed a particular mag- 



netic mechanism for generation of astrophysical gets, a new 
term, "magnetic tower" is gradually gaining popularity too. 
In many computer simulations the development of column- 
like outflows have been observed and these outflows have 
been pronounced magnetic towers, without a detailed analy- 
sis. In fact, coUimated winds and magnetic towers have many 
similar properties like the magnetic braking of the central 
rotator, the amplification of azimuthal magnetic field, the 
axial compression due to hoop stress etc. A magnetic tower 
occupies a finite volume which gradually increases in time 
whereas winds are usually associated with steady-state solu- 
tions that extend to infinity. However, once a wind becomes 
super-fast its interaction with the external medium does not 
affect the upstream solution and thus this distinction is su- 
perficial. A magnetic tower is prevented from lateral expan- 
sion by external pressure but some kind of external force 
is needed to confine any type of magnetic fiow within a fi- 
nite volume. The magnetic field lines of magnetic tower have 
both foot-point anchored to the rotator but the same can be 
said about the wind from a star with, for example, a dipo- 



lar magnetic field; the so-called open magnetic field lines of 
stellar winds close far in the wind zone. 

From our view the main distinction between a mag- 
netic tower and a magnetically driven wind is in the way 
the azimuthal magnetic field is generated in these flows. In 
the magnetic tower model the azimuthal field is generated 
solely due the differential rotation at its base (originally an 
accretion disk) which causes winding of closed magnetic field 
lines whose foot points rotate with different angular veloci- 
ties. Thus, both these foot points are assumed to be in com- 
munication with each other by means of torsional Alfven 
waves; the faster rotating foot-point creates stronger wave 
and wins in the sense that its rotation determines the sigh 
of the generated azimuthal field. This also implies that the 
magnetic tower is a quasi-magnetostatic configuration - the 
vertical growth of the tower can be considered as a sequence 
of magnetostatic models parametrized by time. 

In contrast, wind solutions are characterized by the ap- 
pearance of Alfven horizon; once the Alfven waves emitted 
by the rotator cross this surface they can no longer come 
back. Then the closed field lines of rotator's magnetic field 
may form two groups: (1) those that close within the hori- 
zon and (2) those that close outside of it. If a magnetic field 
line belongs to the first group then its foot points can com- 
municate by means of Alfven waves. In this case it is still 
important whether the rotation at the base is differential 
or not. If, however, it belongs to the second group then the 
communication is broken and for this reason the field line 
can cosidered as opened. For example, the magnetosphere 
of a solidly rotating star with dipolar magnetic field has two 
zones: the so-called "dead zone" where the magnetospheric 
plasma corotates with the star and the magnetic field is 
purely poloidal, and the "wind zone" where the magnetic 
field lines are "opened" and the magnetic field has an az- 
imuthal component. The magnitude of this component de- 
pends on the angular velocity of only one foot point, the 
upstream one, and it is created solely by the Alfven waves 
emitted from this point. One may say that the winding of 
magnetic field in the wind zone is the feedback reaction of 
plasma on the magnetic field that tries to spin it up, it is 
there due to the finite inertia uploaded on the magnetic field 
lines in the wind. There is also finite inertia in the dead zone 
but in the absence of outfiow there is also plenty of time to 
enforce full corotation. 

Formally, the magnetic tower solutions can be con- 
structed for a rather arbitrary distributions of external pres- 



sure (Lynden-Bell 20031. However, for real fiows with finite 



Alfven speed the lateral expansion due to declining external 
pressure can be a factor that works in favor of wind solu- 
tions. Indeed, the non-relativistic Alfven speed Va oc Bp/yfp. 
If zu is the cylindrical radius of a magnetic surface then to 
zero-order approximation we have Bp oc zu~'^ and, assum- 
ing constant mass loading of magnetic field lines, p oc 
This gives Va oc and if the speed of longitudinal ex- 

pansion does not decrease as fast as this then at some point 
the outfiow will become super-Alfvenic. This would signify 
a transition to the wind regime. 

Relativity further complicates the matter as it puts an 
upper limit on the location of Alfven surface - it must reside 



inside the light cylinder (Okamoto 19781. This seems to be 



related to the fact that in relativistic MHD there is inertial 
mass associated with the magnetic (or rather electromag- 
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netic) field itself. In fact, in the limit of Magnetodynamics, 
that is relativistic MHD with vanishingly small particle in- 
ertia (e.g. Komissarov 2002, Komissarov et al.2007), all of 
the inertia is due to electromagnetic field and these surfaces 
coincide - this is why the light cylinder plays such an im- 
portant role in the theory of pulsar winds. In this respect we 
wish to comment on the "pulsar-in-cavity" problem of Uz- 



densky and MacFadyen( 2006| . They consider a magnetized 
neutron star placed in the center of an empty cavity inside 
a collapsing supernova progenitor. The radius of this cavity 
is assumed to exceed the radius of star's light cylinder. By 
making the cavity so large they effectively create conditions 
that allow the development of a pulsar wind with the kind of 
magnetic winding that is specific only to winds (notice that 
the star rotates as a solid body). Rather unfortunately, how- 
ever, they still call the outfiow that may eventually develope 
in this problem a magnetic tower. 

Summarizing, we wish to stress that both in analytical 
and numerical modeling of magnetically driven onflows it 
is very important to determine whether the Alfven horizon 
appears or not. When the magnetic tower approximation is 
used this gives a key test of its validity. When the results of 
computer simulations are analysed this helps to understand 
the physical nature of the outfiow. 



4.4 The nature of magnetar-driven jets 

As the thermally driven explosion (caused by the energy 
deposition in the radiation bubble) develops in our simula- 
tions a rarefaction is created in the center and the magne- 
tar magnetosphere begins to expand, predominantly in the 
equatorial direction. During this initial expansion an Alfven 
horizon appears inside the magnetosphere and a magnetar 
wind begins to develop which injects azimuthal magnetic 
flux and energy into the surrounding space. Because the ex- 
pansion is slow compared to the fast magnetosonic speed 
the fast waves quickly establish approximate magnetostatic 
equilibrium inside the "magnetic cavity" and the hoop stress 
of the azimuthal magnetic field ensures that this equilib- 
rium is characterised by an axial compression. As a result 
the normal stress on the cavity surface develops a maximum 
at the poles, the cavity begins to expand predominantly in 
the polar direction and the polar jets are beginning to form 
(see fig[6|. The jets remain sub- fast up to the very end of 
the simulations and the termination shock never develops. 
However, the jet total pressure is sufficiently high to drive 
bow shocks into the surrounding. In fact the structure of 
onflow is quite similar to that anticipated in Uzdensky and 
MacFadyen pOOe] ). There are however some qualitative dif- 
ferences. First of all, the jets develop before the cavity ex- 
pands beyond the light cylinder. This is not very suprizing 
because the wind magnetization is not ultrarelativistic and 
the Alfven surface can be well inside the light cylinder. Sec- 
ondly, there is no highly relativistic plasma rotation in the 
equatorial plane, even when the cavity expands well beyond 
the light cylinder. 

Eventually, as the central cavity becomes sufficiently 
large, one would expect the magnetar wind to become super- 
fast and the wind termination shock to appear, as it is as- 



cated by the magnetar cooling which leads to higher mag- 
netization of the wind and moves its fast surface further 
away[^ In any case our results suggest that by the time of 
transition to super-fast regime the magnetar-driven jets will 
already be well established and thus the initial setup used 



in (Bucciantini et al. 2007a Bucciantini et al. 2007b I, where 



the central cavity is assumed to be large and spherical, may 
be somewhat unrealistic. 



4.5 Magnetar as a cenral engine of GRB jets 

Our results show that magnetically-driven outflows from 
millisecond magnetars can become highly anisotropic jets at 
the very early stages of supernova explosions. Assuming that 
their propagation speed does not decrease in time, these jets 
would traverse the progenitor star of radius ~ 2 x lO^^cm 
in about 4 seconds. After this the spin-down energy of PNS 
will be carried away by the jets into the surrounding space 
and the rate of energy transfer to the supernova ejecta will 
drop. The total power of long duration GRB-jets is not well 
known but the observations of radio afterglows which are no 
longer subject of rel ativistic beaming su ggest that it can be 
high, ~ 5 X lO^^erg (Berger et al. 20041, well in agreement 



with our results. Since the combined total power of both 
jets in our simulations is ~ 3 x lO^^erg/s, prior to the break 
out the energy transferred to the supernova ejecta will be 
only < lO^^erg, noticeably less than what is usually derived 
for hypernovae. Thus, on one hand our results provide the 
first self-consistent numerical model of a central engine capa- 
ble of producing, on a long-term basis, coUimated jets with 
sufficient energy to explain GRBs and their afterglows. On 
the other hand, the relatively small amount of energy trans- 
ferred the supernova ejecta could be a problem due to the 
connection between high velocity SNe and long GRBs. Al- 
though the connection between SNe and GRBs is supported 
by strong circumstantial evidence, the direct evidence is still 
scarce with only handful of SNe identified with GRBs. It is 
quite possible that the current paradigm connecting GRBs 
with HNe is simply a reflection of the observational bias to- 
wards more powerful events ( [Woosley fc Bloom 2006| |Della| 
Valle2006[ ). SN 2002LT/GRB 021 211 may the case o f a GRB 
produced by a standard SN Ib/c ( |Defla Valle 2006 1. 

The relatively high propagation speed of jets in our 
simulations may in part be attributed to the condition of 
axisymmetry. The hydro-simulations of 2D slab jets (e.g. 
Komissarov & Falle 2003) and 3D round jets ( Norman 1996 1 



sumed in ( Bucciantini et al. 2007a Bucciantini et al. 2007b I 



However, additional studies are required to determine when 
exactly this will occur. This issue may be further compli- 



show signiflcantly lower propagation speeds of jets without 
enforced mirror or axial symmetry as the kink modes of 
current-driven or pressure-driven instabilities result in redis- 
tribution of the pressure force over larger area at the jet head 
(cf. Aloy et al.,1999). On the other hand, as the jets enter 
progressively less dense outer layers of collapsing star their 
propagation speed may actually increase. Slower propaga- 
tion speed would increase the amount of energy transferred 
to the supernova ejecta and reduce the energy available to 
produce GRB. 

Although our jets are not as fast as required by the ob- 
servations of GRBs the gradual cooling of PNS would lead 
to lower mass loading and hence higher terminal speed of 



For inertia free wind this surface moves to infinity 
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the jets later on. Indeed, according to the computations of 



Pons et al.(1999l the neutrino luminosity of PNS decreases 
by a factor of 10 already during the first 10s and then en- 
ters the phase of rapid exponential decline, largely indepen- 
dent on the details of initial models and EOS. Metzger et 
al.(2007l argue that the magnetization a = 100 — 1000 can 



be reached when the the PNS still has enough rotational 
energy to power a GRB. This corresponds to the asymp- 
totic Lorentz factor in the range of Foo = a/2 = 50 — 500 
(e.g. Vlahakis 2004). However, the magnetic acceleration of 
relativistic fiows is a rather slow process and we do not ex- 
pect the termination speeds to be reached within the stellar 
boundaries. Further studies are needed to clarify the issue 
of magnetic acceleration. 

Since the cooling time of PNS is comparable with the 
break out time, and the time-scale of long duration GRBs, 
one would expect that that at the time of break out the 
mass-loading of the PNS magnetosphere will still be quite 
high and thus the flow speed of the jets will still be relatively 
low, F ~ 1 — 10. Only later the Lorentz factor will gradually 
increase and eventually reach the ultra-relativistic values, 
F ~ 200, inferred for the GRB jets. The implications of 
such a non-uniform jet structure, with Lorentz factor grad- 
ually decreasing with distance from the the star, for both 
the prompt and the afterglow emission of GRBs remain to 
be investigated. Here we point out only few obvious points. 
First of all one may expect some similarities with the model 



of structured jet (Meszaros et al. 1998 


Dai & Gou 2001 


Panaitescu 2005| |Rossi et al. 2002| |Wei & Jin 2003| |Kumar| 


& Granot 2003 


1, whose Lorentz factor varies with polar an- 



gle, and the model of two-component jet ( 


Meszaros & Rees 


2001| Vlahakis et.al. 2003| Zhang et al. 2003| Ramirez-Ruiz 


et al. 2002 Mizuta et al. 2006||Granot et al. 2006| Jin et al. 


2006| Morsony et al. 2006| Willingale et al. 2007 


1. Secondly, 



since lower Lorentz factor results in weaker Doppler effect 
one would expect softer emission from the parts of the jet 
emitted earlier and located further away from the star. This 
could be the origin of the prompt optical and X-ray emission 
and early afterglows. Curiously enough, the model suggests 
that the zero-point of gamma-ray bursts lags behind the zero 
points of prompt X-ray and optical bursts. Finally, while 
the collision between the slower earlier jet and the ISM or 
the progenitor wind would still produce the strong forward 
shock usually associated with afterglows, one might expect 
a secondary strong forward shock where the faster late jet 
collides with the slower early jet. This secondary shock will 
propagate faster than the primary one and will eventually 
catch up with it. Its emission, which will be harder and 
beamed more strongly than that of the primary shock, may 
lead to distinctive features in the light curves of afterglows. 

It is tempting to consider the global oscillation of 
magnetar's magnetosphere and the related non-stationary 
plasma ejection as the origin of "inner shocks" invoked in 



models of prompt GRB emission (Piran 2005b I. However, 
it is important to check if the reconnection rate is deter- 
mined mainly by the global dynamics and not by the resistiv- 
ity model (which was purely numerical in our simulations). 
Moreover, it remains to be seen if such oscillations can per- 
sist at later times (t > few seconds) when the the mass- 
loading of the magnetosphere drops and the inner "cavity" 
significantly increases in size. Later, when the jet becomes 
super-fast, internal shocks can be produced via interaction 



with inhomogeneous time-dependent cocoon and instabili- 
ties. The inhomogeneous structure of the slow jet may also 
lead to variable emission from the secondary forward shock. 



5 CONCLUSIONS 

The results of our study show that when a core collapse of a 
massive star results in a birth of a millisecond proto neutron 
star with super-strong magnetic field (proto-magnetar) , this 
could have a spectacular effect on the supernova explosion. 
During the very early stage of the explosion the magnetar 
can produce highly collimated jets capable of puncturing 
the collapsing star in a matter of seconds. The spin-down 
time of the proto-magnetar is an order of magnitude longer 
thus suggesting that a large fraction of its rotational energy, 
E ~ lO^^erg, will be carried by the jets into the surrounding 
space. This supports the idea that at least some long dura- 
tion Gamma Ray Bursts can have millisecond magnetars as 
their central engines. 

The magnetic outflow is best described as a sub-fast 
super-Alfvenic collimated wind. The super- Alfvenic nature 
of the flow explains the generation of azimuthal magnetic 
field which soon begins to dominate the fiow dynamics. The 
collimation is a combined effect of the inertial confinement 
by the stellar material and the hoop stress of the azimuthal 
field. The outfiow is not a magnetic tower. 

It remains to be seen as to how soon the wind in the cen- 
tral cavity becomes super-fast and a proto-PWN is formed 
inside the collapsing star. The magnetic acceleration of the 
jets inside the channels bored through the star and outside 
of the star is also an important subject for future study. 

It is very likely that the GRB jets produced in this 
fashion first emerge as only moderately relativistic fiows 
and only later as the magnetar cools they become ultra- 
relativistic. The effect of this on both the prompt and the 
early afterglow emission needs further investigation 

A non-magnetic, e.g. the delayed neutrino-driven ex- 
plosion of power comparable to that of normal core-collapse 
supernovae is needed to turn-on the magnetic mechanism. 
Otherwise the PNS magnetosphere is unable to expand end 
develop a super-Alfvenic wind. The idea that a failed su- 
pernova explosion can be revived by a millisecond magnetar 
does not seem to work. 
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